Parallel pumping of electrons 
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We present simultaneous operation of ten single-electron turnstiles leading to one order of mag- 
nitude increase in current level up to 100 pA. Our analysis of device uniformity and background 
charge stability implies that the parallelization can be made without compromising the strict re- 
quirements of accuracy and current level set by quantum metrology. In addition, we discuss how 
offset charge instability limits the integration scale of single-electron turnstiles. 



Realization of a standard for electrical current based 
on discreteness of the electron charge e is one of the ma- 
jor goals of modern metrology. The theoretical basis for 
obtaining current I = nef when n electrons are sequen- 
tially transferred at frequency / has been well known for 
more than two decades [Jllli- With multijunction de- 
vices it has been possible to demonstrate pumping with 
relative accuracy of 10~ 8 up to picoampere level (4{. Al- 
though the accuracy is more than an order of magnitude 
better than the present definition of ampere, the out- 
put current level is, however, too small for applications 
apart from the capacitance standard [5j . In order to cre- 
ate large enough current in this fashion with the desired 
accuracy, parallelization of multiple pumps is inevitable. 
In this letter, we demonstrate parallel electron pumps 
with quantized current plateaus. The parallelization is 
done up to ten devices leading to a current level exceed- 
ing 100 pA. This is already enough for the closure of 
the so-called quantum metrological triangle [y, |2] which 
would then justify the current standard based on single 
electron transport. 

The idea in quantum metrological triangle is to probe 
the consistency of the current from an electron pump 
against two other quantum phenomena, resistance from 
Quantum Hall effect and voltage from the AC Joseph- 
son effect. This verification would yield a consistency 
check for two fundamental physical constants, the charge 
of electron e and Planck's constant h, and enable one to 
define the Si-units of electrical quantities directly from 
quantum mechanics. To obtain higher current levels, 
various approaches have been studied @, 0, [l(J U, 12, 
13, LLJ], such as surface acoustic waves, superconducting 
devices and semiconductor quantum dots but still, the 
present accuracy of these devices is limited. Two separate 
semiconductor quantum dot devices have recently been 
operated in parallel [la ]. The hybrid turnstile [HI, [13] 
used in this work holds the promise of achieving ex- 



tremely low pumping errors 18|, similar to the multi- 
junction circuits. In addition, due to the simplicity, the 
turnstiles can be scaled up to higher integration levels for 
parallel operation as they require only one tuning signal 



per device. 

The scheme of parallel turnstiles is shown in a scan- 
ning electron micrograph in Fig. [T] (a). The samples 
used in this letter were fabricated with standard elec- 
tron beam lithography on top of a spin-on-glass insulator 
layer as explained in supporting information. In each of 
the repeated cells there is one individual device. It is a 
single-electron transistor (SET) where the tunnel junc- 
tions are formed by an overlap between superconducting 
leads and a normal metal island. With respect to paral- 
lelization, these devices require one independent DC gate 
voltage V 9j i per device for compensating the inevitable 
offset charges. The other signals, bias voltage Vb for set- 
ting the preferred tunneling direction and the RF gate 
voltage Vrf used for pumping can be common for all the 
devices. 

The operation of a hybrid turnstile can be under- 
stood by considering the energy thresholds for single elec- 
tron tunneling. The thresholds are determined by the 
externally controllable bias voltage Vb and gate offset 
n g = (C gii V g> i + C RF ,iV RF )/e, where C g ^ and C RF ,i are 
the coupling capacitances from the DC and RF gates to 
the island respectively. In Fig. Q] (b), we present the 
measured DC current and the thresholds of one of the 
devices on the Vb - n g plane. The stable regions for the 
charge states n = and n = 1 (red and black boxes 
respectively) overlap, and ideally no DC current flows 
even at a finite bias voltage up to V> = 2A/e. Here A 
is the energy gap of the superconductor and n the ex- 
cess number of electrons in the island. The broad stabil- 
ity regions enable one to pump electrons sequentially by 
moving the gate offset n g along the horizontal pumping 
trajectory shown as the solid blue line. To obtain high 
accuracy, the turnstile should spend enough time beyond 
the thresholds shown as solid lines but should not cross 
the thresholds of backtunncling shown as dashed lines. 
All devices should cross the forward tunneling thresholds 
in concert while avoiding the backward tunneling. 

To test such uniformity in a parallel set-up, four turn- 
stiles were connected to a common drain while the source 
sides were left separate for individual characterization. 
Afterwards they were connected together to a common 
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FIG. 1: Parallelization scheme and the operation principle of a single turnstile, (a) Scanning electron micrograph of parallel 
turnstiles. The figures are coloured to increase clarity. Normal metal islands (red) are placed on top of a common radio 
frequency (RF) gate (yellow). This gate is insulated from the islands with a SiO x layer. Blue regions denote the superconducting 
aluminium wires. The tunnel junctions are formed under the islands by oxidising aluminium. One gate line with DC voltage 
V g ,i is needed for each of the devices while the bias voltage Vb and the RF signal Vrf can be common to all turnstiles, 
(b) Stability diagram of the turnstile D (See Table [TJ) and the tunneling thresholds for electron pumping. For tunneling 
into the island the requirement is —2E c (n + 1/2 — n 3 ) ± eVb/2 > A and for tunneling out from the island it is given by 
2E c (n — 1/2 — n 9 ) ± eVb/2 > A. The upper (lower) sign corresponds to the junction which lies on the positive (negative) side 
of the bias. E c — e 2 /2C is the charging energy of the island with total capacitance C and n the number of excess electrons in 
the island. n g — (C g ,iV gt i + Crf^Vrf) /e is a normalized gate induced offset charge used for controlling the energy thresholds. 
The green region inside solid black lines is stable for n = while the green region inside solid red lines is stable for n = 1. Solid 
lines are thresholds for wanted transitions during pumping while dashed lines correspond to backtunneling in wrong direction. 
Solid blue curve shows ideal pumping curve with positive bias voltage at gate open, (c) Current-voltage characteristics of the 
turnstile D. The gate offset charge is swept back and forth so that the envelopes correspond to gate being open or closed. 
Simulations used to extract device parameters from these extremum cases are shown by black solid lines. 



bias voltage source to demostrate parallel pumping. The 
DC current-voltage characteristics of one turnstile is 
shown in Fig. [T] (c) where the gate voltage is swept 
back and forth so that we obtain both extreme cases of 
gate open (n g = 0.5) and gate closed {n g = 0.0). Also 
the current-voltage simulations for both of these cases 
are shown. These curves are calculated with sequential 
tunneling approximation and used to extract the device 
parameters which are listed in Table [T] From pumping 
measurements, the rising edge to the first plateau was 
determined as the RF gate voltage V r , m for which the 
current is half of the value at plateau. The measured vari- 
ation of V rym between the devices was ±7 %. According 
to numerical simulations, this narrows the metrologically 
fiat part of the plateau by about 10 %. 

In addition to the crucial parameters determining the 



TABLE I: The parameters of the four turnstiles A-D. Rt 
and Crf are estimated from the measurement data with un- 
certainty of 1 %. A and E c /A are fitted with the help of 
numerical simulations to within 2 % precision. 



Rt (k») A (/iV) ffc/A Crf (aF) V r , m (mV) 



A 


490 


213 


1.03 


25.3 


3.82 


B 


580 


214 


1.10 


23.5 


3.60 


C 


610 


214 


1.10 


24.7 


3.83 


D 


742 


215 


1.16 


23.4 


4.12 



thresholds, also individual tunneling resistances Rt of 
the turnstiles are obtained from simulations. This pa- 
rameter together with the total capacitance C and su- 
perconducting energy gap A determine the maximum 
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FIG. 2: Parallel pumping of four turnstiles at / = 10 MHz. From (a) to (d), one to four devices are tuned between gate open 
and gate closed states respectively while the rest of the turnstiles are kept at gate open. This yields plateaus where 4 or 4 ± N 
electrons are pumped in each cycle where N is the number of tuned devices. The case N = 4 shows pumping curves similar 
to those of a single turnstile but with four times higher current. The slope on the plateau where all devices are pumping one 
electron is 810 Gtt per turnstile with respect to RF amplitude. The slopes for single turnstiles were similar up to the accuracy 
of this measurement limited by the current drift on the 10 fA level. 



operation frequency of a turnstile. However, for paral- 
lellization, there are no constraints on the similarity of 
the tunneling resistances. Yet, the largest of them deter- 
mines the maximum operational frequency of the system. 
For an aluminium based device, the maximum current is 
limited to somewhat above 10 pA when metrologically 
accurate operation is required fljl . Il9j |. Therefore the 
parallel operation is inevitable to obtain higher current 
levels while simultaneously preserving high pumping ac- 
curacy. 

After the characterization of individual turnstiles, they 
were connected in parallel and pumping curves presented 
in Fig. [5] were measured. Here we have changed the gate 
states of one (Fig. [2] (a)), two (b), three (c) and four 
(d) devices simultaneously while keeping others at gate 
open. Thus we obtain current plateaus where zero to 
eight electrons are transported within one cycle. This 
measurement demonstrates that we can fully control the 
DC gate states of each device. The gate state of each 



individual device i was extracted from the total current 
through all of the devices by sweeping gate voltage V g ^. 
The gate open states correspond to maximum values of 
current. Cross-coupling between the gates was less than 
3 % and hence only one iteration round after a rough 
setting of the gates was needed to get gate states correctly 
within 1 %. 

Next, to demonstrate the reproducibility and robust- 
ness, ten turnstiles were operated similarly by a single 
RF drive. Two chips were used from different batches 
with six turnstiles on one chip and four on the other. All 
ten devices were bonded to one common bias line and 
hence no preliminary characterization of individual de- 
vices was made. The results for different bias voltages 
are shown in Fig. [3] This setup yields 104.1 pA at the 
first plateau with a pumping frequency of 65 MHz, which 
demonstrates a current level large enough for closing the 
quantum metrological triangle [7J. In the present exper- 
iment, the number of parallel devices was limited by the 
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FIG. 3: Ten turnstiles working at / = 65 MHz. All turnstiles 
are set either to gate open or gate closed state. The current at 
the first plateau is approximately 104.1 pA. The insensitivity 
to the bias voltage is also shown. The numbers in the panel 
indicate bias voltage values. 



is somewhat similar to previously reported values 2^, 21 1 . 
We note that even better performance with no drifts has 
been observed for metallic single-electron devices previ- 
ously. 2^, 23|. Such improvement would allow further 



increase in the number of parallel devices. 

The main result of the present work is the controlled 
operation of parallel electron pumps. The stability of 
offset charges was studied and it is shown to allow more 
than ten parallel devices to be operated without signif- 
icantly compromising the accuracy. These devices are 
prominent to fulfill the strict accuracy requirements for 
the closure of the quantum metrological triangle, and as 
the outcome of this work, we show that the obtained 
current level achieves the generally accepted goal. The 
flatness of the plateaus ispreserved in parallel operation. 



Moreover, Refs. [17|, \2M an d our recent unpublished 



work suggest that the sub-gap leakage, which is the re- 
maining error source, can be significantly decreased by 
proper design of the electromagnetic environment. 
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number of DC lines available. 

In more general terms, the number of devices which 
can be operated in parallel simultaneously is determined 
by the offset charge stability. The strategy of a pump- 
ing experiment is to first set the gate of each turnstile, 
then perform the pumping measurement and afterwards 
check the offset charges again. If they are not within the 
limits, one discards the data. To study the stability of 
the ten devices, the DC gate modulation was measured 
as a function of time for each of the turnstiles simul- 
taneously. The time for one cycle was chosen to be 10 
minutes which was equal to the time required to measure 
the data of one panel in Fig. [21 The gate stability for a 
typical turnstile is shown in Fig. [4] (a). From the mea- 
sured data, histograms of the offset charge changes were 
determined for each of the devices as shown in Fig. [4](b). 
Moreover, in the inset a histogram for the corresponding 
maximum change of the ten turnstiles is presented. From 
this we obtain 73 % probability to get valid data with this 
measurement setup as described in caption. Addition- 
ally, we can estimate the maximum number of turnstiles 
operable in parallel to be 17 in the present case, which 
would yield efficiency of 50%. By making the measure- 
ment period smaller one could increase the number of 
devices as they have less time to get offset. We estimate 
that the measurement period can be decreased by one or 
two orders of magnitude. This will allow one to increase 
the number of parallel devices accordingly. Also differ- 
ent materials or fabrication methods can provide smaller 
drifts and hence allow larger integration scale. In our 
devices, typical spectral density of charge noise followed 
relation S q (f) = a 2 / f 2 at the observed frequency range 
/ = 1 /iHz — 1 niHz with a = 10 -6 e\/Hz. The magnitude 
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SUPPLEMENTARY INFORMATION, METHODS 
USED IN SAMPLE FABRICATION AND 
MEASUREMENTS 

Our circuits comprising parallel turnstiles were fabri- 
cated at NEC Nano Electronics Research Laboratories, 
Tsukuba Japan, on thermally oxidized Si chips using 
one photolithography and three electron-beam lithogra- 
phy steps. The chip size is 3.6 mm x 3.6 mm. First, 
Ti/Au (5 nm/95 nm) leads and contact pads were made 
using a standard liftoff photolithography process. The 
pads are located at the edges of the chip with the leads 
stretching to the chip center leaving a square-shaped 80 
/im x 80 /im area free for fine structures. Then, in this 
area, a 3 nm / 20 nm thick Ti/Au RF gate was deposited 
through a conventional soft mask formed in bi-layer re- 
sist by electron beam lithography. The RF gate was con- 
nected to one of the leads and covered in the next step by 
a patterned spin-on glass layer to isolate it from the turn- 
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FIG. 4: Stability of offset charges, (a) Gate modulation of one of the turnstiles as a function of time. The light/dark areas 
correspond to maximum/minimum current. The red dots show one of the gate open states which is changing due to variations 
in the offset charges. This data is obtained as one set within sequential sweeping of the gate charges of the ten turnstiles. These 
ten sweeps are then repeated every ten minutes, (b) A histogram of the changes in the offset charge between each sweeping set. 
The different colours denote different individual devices. In the inset, a histogram of the largest change in each set is shown. 
Total number of sweeps per turnstile is 410. Although for individual devices the changes are peaked at zero, for ten turnstiles 
it is more likely that at least one of the devices has a few percent change. To obtain an estimate for the maximum number 
of devices that can operate in parallel, we assume that the offset charge changes are independent as the average correlation 
coefficient between the devices was 0.12. This gives a lower limit for the probability of having N turnstiles in correct gate 
states as pn = Pi , where pi is the probability of having one turnstile in a correct state. The flatness of the theoretical current 
plateaus is such that with the observed 10 % variation in tunneling thresholds we can still tolerate 5 % changes in gate offsets. 
With the measured data, this will lead to pi — 0.96 and pio = 0.71. This is consistent with the value pio = 0.73 obtained from 
the data in the inset of (b). Therefore the requirement for efficiency of pjv > 0.5 will limit the number of parallel turnstiles to 
N — 17 according to the presented data. 



stiles made above the RF gate. Finally, SINIS type turn- 
stiles were fabricated using two-angle deposition though 
a suspended mask created in a Ge layer using tri-layer 
electron-beam process. The turnstile pattern is exposed 
in the top layer of polymethylmetacrylate and then, af- 
ter development, transferred into Ge layer by reactive ion 
etching in CF4 gas. The undercut under Ge mask was 
formed by etching of the bottom copolymer with oxygen 
in the electron-cyclotron-resonance machine. Deposition 
of the turnstile leads (Al) and islands (Au/Pd) was done 
immediately before the low temperature measurements 
in an e-gun evaporator with a oxidation step in between. 



Measurements were performed at Low Temperature 
Laboratory at Helsinki University of Technology in a 
3 He/ 4 He dilution refrigerator at temperatures below 
0.1 K. Standard DC and RF voltage sources were used 
with resistive dividers and attenuators to set the op- 
eration voltages. Current was measured with a room- 
temperature low-noise current amplifier. 
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